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The electroweak 0{a) radiative corrections to the decay of the heavy CP-even MSSM 
Higgs boson to weak gauge bosons are presented. Due to the suppression of the tree-level 
H^WW / Z Z coupling, the electroweak contributions to the partial decay width are signif- 
icant. Although the effective Born decay width can be rather small for certain M^o values 
(especially for large tan/3), the corrected partial widths at different values of M^o are of 
Oh| comparable size. 

1 Introduction 

One of the major tasks of the Large Hadron Collider (LHC) at CERN is to explore the mechanism 
■ responsible for the electroweak symmetry breaking (EWSB). In the standard model (SM), the 



OO 



o 



X 



EWSB is realized by the Higgs mechanism, which predicts the existence of a physical scalar boson, 



0^ ' the Higgs boson. The mass of the Higgs boson cannot be predicted by the SM. Experimental 



searches at LEP [Ij and Tevatron [2] have excluded a Higgs mass Mh < 114.4 GeV and 158 GeV < 
Mh < 173 GeV, both at 95% confidence level (C.L.). In addition, electroweak precision analysis 
favors a relatively light SM Higgs boson with a mass Mh < 158 GeV (or Mh < 185 GeV if the 
LEP2 direct search limit is included) [3j. 

In the minimal super symmetric standard model (MSSM), which is the most economic exten- 
sion of the SM that incorporates supersymmetry, two complex Higgs doublets are required for 
consistency. After the EWSB, three of the eight degrees of freedom are absorbed by the weak 
gauge bosons and become their longitudinal components, leading to five physical Higgs bosons, 
/i° and H° (CP-even), A° (CP-odd) and (charged). At tree-level the Higgs sector of MSSM 
can be described by two parameters, which can be chosen as M^o and tan (3, where M^o is the 
mass of the CP-odd Higgs boson and tan /3 is the ratio of the vacuum expectation values of the 
two Higgs doublets. The tree- level mass of the light CP-even Higgs boson is bounded from above 
by the mass of Z boson as a consequence of supersymmetry. Dependence of the Higgs sector on 
parameters of other sectors enters via radiative corrections. By including radiative corrections 
(up to two-loop order), the upper mass bound of the light CP-even Higgs boson is shifted to 
~ 135 GeV [IIIIIIH]- 
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It is interesting to investigate the behavior of the MSSM Higgs sector in the decoupUng Umit 
Mj^o ^ Mz [9]. In this hmit, the heavy CP-even Higgs boson and the charged Higgs bosons 
are nearly degenerate in mass with A^, and the coupling of the light CP-even Higgs boson to 
SM fermions and gauge bosons resembles the corresponding coupling for the SM Higgs boson. If 
only one light Higgs boson is discovered, it might be compatible both with the SM and with the 
MSSM. If there are extra heavy Higgs bosons observed, they would indicate the incompleteness of 
the SM. Investigating the decay properties of such heavy Higgs bosons can then help disentangle 
supersymmetry from other potential extensions of the SM with an extended Higgs sector. 

In the MSSM, the tree-level coupling of the heavy CP-even Higgs boson to gauge bosons is sup- 
pressed by a factor of cos(/3 — a), compared to the corresponding coupling for the SM Higgs boson, 
where a is the angle that diagonalizes the CP-even Higgs sector at tree-level. This suppression 
can be rather strong for large values of My^o, hence it is of particular interest to investigate the 
impact of radiative corrections in such situations. 

In this paper we consider the electroweak radiative corrections to the decay of the heavy CP-even 
Higgs boson to weak gauge bosons. Electroweak corrections can induce important modifications 
to the tree- level H^WW/H^ZZ coupling. One potential source of large corrections is the contri- 
bution of loops involving fermions and sfermions, and in particular of loops involving the third 
generation fermions and sfermions, since they contain potentially large Yukawa couplings. The 
Higgs propagator corrections can give rise to significant contributions as well. In this work we 
concentrate on CP-conserving MSSM with real parameters, therefore the heavy CP-even Higgs 
boson can only mix with the light CP-even one. 

The outline of this paper is as follows. In section[2]we discuss the structure of the H^WW/H^ZZ 
vertex correction. In section [3] the decay amplitudes of to WW/ZZ are given. Section U] is 
devoted for numerical discussions. We draw our conclusions in section [5l 



2 Corrections to the H^WW/H^ZZ vertices 

2.1 Correction to the H^WW vertex 

The tree-level H^WW coupling is given by 

= ^ cos(/3 - a)g^^- ^ Vh^ cos(/3 - a)^^^ , (1) 

where V^o denotes the coupling of the SM Higgs boson to W bosons. When the mass of the CP- 
odd Higgs boson M^o becomes large, the angle /? — a — )■ 7r/2 and the factor cos(/3 — a) approaches 
as 0(M|| sin4/3|/M^o)j thus the tree- level coupling is strongly suppressed. 

The one-loop corrected H^WW vertex possesses the following structure [11] 

= y^o,o + ^i^O (^^2 ^2 + ^^3 ^3 + C^2 ^3 + ^^3 ^2 + Eg^"" - iFe^^f"' k2pha) , (2) 

where k2, denote the momenta of the two gauge bosons, et^'^P'^ is totally antisymmetric with 
g0i23 _ j£ gauge bosons are on-shell, then only D and E terms will contribute. As 
mentioned in the introduction, the loop contributions from the fermionic and sfermionic sector, 
and in particular from the third generation fermions and sfermions are expected to be sizable. In 
Fig. [1] we show as examples the Feynman diagrams involving fermion and sfermion loops. The 
packages FeynArts [13] • FormCalc [TB] and LoopTools [TO] are used throughout the computation. 
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Figure 1: One-loop vertex diagrams and counter term diagram for H^WW . 
The contribution of the vertex counter term diagram in Fig. [T]can be written as (see e.g. |20j ) 
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6sw , sin(/3 - a) 2 or, « ,1 r-z x 

+ + T3 r(-cos /3 5tan/3 + -5ZhOHo) 

sw cos(p — a) 2 



(3) 



where dZ^ is the charge renormahzation constant, 6Z\\r is the field renormalization constant for the 
W boson, and 5My^ the mass counter term for the W boson. 5sw is the renormalization constant 
for the weak mixing angle. SZj^ouo and SZ jjojjo are the Higgs field renormalization constants in 
the mass eigenstate basis. The counter term for tan /? is introduced via tan /3 — )• tan (3 + 6 tan /3. 
The vertex counter term is proportional to the tree-level coupling, and thus contributes to the 
^-term in Eq. 



To determine the Higgs field renormalization constants, we choose the DR scheme [2T|. In this 
scheme, the Higgs field renormalization constants are given by 

5Z;,o;,o = -[ReS;o;,oKo)]div, 

5ZjjOfjo = — [ReSj 



5Zj-fijj[0 = 5Zjj[0 



sin a cos a 



cos 2a 



{SZ^OhO - 5 Z HQ ho) , 



(4) 



where m/jO, nijjo denote the tree-level masses of the two CP-even Higgs bosons, "div" means we 
only take into account the divergent parts of the renormalization constants. The counter term 
6 tan f3 is fixed by 



5 tan /3 f 5 tan /3 



tan/3 



tan /3 



DR 



1 



2 cos 2a 



{dZ^OhO - SZhoho) ■ 



(5) 



For the renormalization scale, we choose //^^ = mt- The remaining counter terms in Eq. ([3]) are 
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fixed in the on-shell scheme as follows |22j 



sw 
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sw 



cw 

-ReS'^(M^^ 
ReE^(M^) , 
ReSi(M|) , 



2 s2 



(6) 



where the prime indicates the derivative and the superscript T denotes the transverse part of the 
corresponding self energy. 

Due to the mixing between the two CP-even Higgs bosons beyond tree-level, when evaluating 
the radiative corrections to the decay WW we need to take into account the h^WW vertex 

correction as well. The tree-level h^WW coupling is simply obtained by replacing cos(/3 — a) in 
Eq. ([T|) with sin(/3 — a), and the corresponding counter term contribution reads 



Vfjo sin(/3 — a)g 



R7 ^l^^^H^ 



^ bsw ^ cos(/3 



sin(/3 



— -(cos^ /? Jtan /3 



(7) 



where the field renormalization constants SZf^OfiO and dZfjoj^o are given in Eq. 



2.2 Correction to the H^ZZ vertex 

The tree-level H^ZZ coupling has the same structure as the tree-level H^WW coupling 



yflU 



eMw 
Cwsw 



cos(/3 - a)gf"' = V^o cos(/3 - a)g'"' , 



(8) 



it also differs from the SM counterpart by a factor of cos(/3 — a) and gets strongly suppressed in 
the decoupling limit M^o » Mz- 

The corrected vertex has the same tensor structure as in Eq. ([2]), with Vfjo replaced by V^q. 
The corresponding counter term contributions are given by 



X7 , X7 , ^ ^^W ^sw .^s'^ sin(/3 - a) 2«., n 

dZe + oZz + - , ^9 \^~T — 1) + TTTTS 7 (—cos p o tan p 
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and 



Kr^^ = V^jo sin(/3 - a)g^^ 



(9) 
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with 6Zz = -ReS'|(M|) 
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3 Decay amplitudes of WW/ZZ 

The decay amplitude can be obtained from the couphng vertex in the previous section by putting 
all external momenta on-shell and multiplying with the polarization vectors for the external 
gauge bosons. Beyond the lowest order, the CP-even Higgs propagator matrix receives important 
radiative corrections, leading to finite wave function normalization factors for the external Higgs 
boson. These wave function normalization factors can be incorporated by using the following 
effective Born amplitude 



hOJ 

= v^-MV(1 + tan(/3 - a)Zj^o^o) , (11) 

where a-nd M^o denote the tree-level decay amplitude for and h^, respectively. For 

illustrating purposes, we will present in the next section both the tree-level results obtained from 
A^^o and the effective tree-level results obtained from -M^q- The wave function normalization 
factors Zffo and Zfjo^o in Eq. ()lip can be determined from the renormalized self energies of Higgs 
bosons as 



Zho 



1 



1 + ReS' o(^^) - Re 



where Mfjo is the physical mass of H^. In this work the physical masses of Higgs bosons and the 
finite wave function normalization factors are computed with the program package FeynHiggs [23] , 
in which the dominant two-loop corrections to the Higgs boson self energies are also taken into 
account. For the decay to W bosons, is given by 

= yROfi'i^^- ' (13) 

where V^o q is defined in Eq. ([T]) and e^^^ are the polarization vectors of the external W bosons. 
The amplitude for the decay of to Z bosons can be obtained analogously. 

At one-loop level, one can also write an effective amplitude 

A^eff = Vz^o{M\jo + ZHOf^oMlo) . (14) 

In the following the results obtained from this effective amplitude will be denoted as effective 
one-loop results. In obtaining such results we also include the square of the one-loop amplitude, 
since the tree-level coupling can be suppressed so that the square of the one-loop amplitude can 
become comparable to the tree-level result. For the decay WW, we compute the complete 

0{a) contributions, and the effective one- loop contributions from the fermionic and sfermionic 
loops, which do not involve infrared divergences. For the decay — )• ZZ, we compute also the 
complete effective one-loop contribution as the complete one-loop amplitude is infrared finite in 
this case. 



4 Numerical discussions 

For the numerical evaluation, we choose the benchmark scenarios suggested in [26], where the two 
parameters that govern the tree-level Higgs sector, M^^o and tan/?, are kept as free parameters. 
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The SM parameters used in the numerical analysis are the same as in ref. [28] , Throughout the 
parameter scan, the experimental mass exclusion limits from direct search of supersymmetric 
particles and the upper bound on the SUSY corrections to the electroweak p parameter [29j have 
been taken into account (except for the tree- level result). In the parameter region 50GeV < 
Mjifi < 1 TeV and 1 < tan /3 < 50, the gluophobic scenario has been ruled out by the bound 
derived from the BR(i? — ?> ^^7) prediction [30], therefore we will not discuss this scenario here. 
The investigated scenarios are listed as follows: 

1. The m^^^ scenario 

The parameters in this scenario are given by 

MsusY = 1 TeV , = 200 GeV , M2 = 200 GeV , 

Xt = 2MsusY , Ab = At = Ar, mg = O.SMgusY , (15) 

where MgusY is the soft SUSY-breaking parameter, // is the supersymmetric Higgs mass param- 
eter, M2 denotes the SU{2) gaugino mass, Xt the mixing parameter of the top squark sector, 
Afj^t^T the trilinear couplings for the third generation squark and slepton and nig the gluino mass. 
This scenario yields a maximal value of the lightest CP-even Higgs boson mass for given M^o 
and tan /3. 

2. The no-mixing scenario 

The only difference of this scenario from the m]!^^'^ scenario is the vanishing mixing in the top 
squark sector and a higher value of MgusY, where the latter is chosen to avoid the exclusion 
bounds from the LEP Higgs searches [HEI]- The parameters in this scenario read 

MsusY = 2 TeV , /i = 200 GeV , M2 = 200 GeV , 

Xt = , Ab = At = Ar , mg = O.SMsusY • (16) 

3. The small-Ogfj scenario 

In this scenario a suppression of the h^bb coupling can occur. The parameters are given by 

MgusY = 800 GeV , /i = 2.5Msusy , M2 = 500 GeV , 

Xt = -1100 GeV , Ab = At = Ar, m-g = 500 GeV . (17) 

To illustrate the numeric impact of radiative corrections, we show in Fig. [2] and [3] the dependence 
on Mj^o of the lowest order and the corrected partial decay widths in the m^^^ scenario for 
tan /? = 5, 30. As can be seen by comparing the tree-level and the effective Born results in 
the figures, the contributions of Higgs propagator corrections are significant, they change the 
dependence of the tree-level result on M^o dramatically. An interesting feature of the effective 
Born results is that they reach a minimum for moderate M^o values (~ 420/500 GeV depending 
on tan /3). This is due to the cancellation of the two parts in the effective Born amplitude Eq. (|lip 
at such values of M^o. As shown in the figures, the tree-level partial decay width for tan/3 = 30 
is much smaller than that for tan/3 = 5, since for relatively large Mj^o, the tree-level coupling 
of to vector boson pair is suppressed by tan /? as well, as explained in the discussions below 
Eq. ([1]). For — > WW, at tan/3 = 5 the contribution of the fermionic and sfermionic sector 
yields the dominant part of the 0{a) corrections, and the leading contribution from the fermionic 
and sfermionic sector is that from the third generation fermions and sfermions. For tan /3 = 30, 
besides the fermionic and sfermionic contribution, the contribution from other sectors to 0{a) 
corrections are also important. In Fig. [2] we also show the effective one-loop results with fermionic 
and sfermionic loop contributions (corrections beyond 0(a) described in SecOare included), from 
which one can see again that the third generation fermionic and sfermionic contribution yields the 
leading contribution of the fermionic and sfermionic sector. For tan /3 = 5, the relative difference 
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between the effective one-loop results and the 0{a) ones is less than 50% in the plotted region. 
For tan/3 = 30, the relative difference between the effective one-loop results and the 0{a) ones 
is less than 50% only for small to moderate My^o values, it can go beyond 100% for large M^o 
values and reaches ~ 200% for M^o = 800 GeV, as a consequence of the strong suppression of 
the tree-level coupling by both tan f3 and Af^o . 

Fig. [3] depicts the partial decay width for — > ZZ. The decay width falls off rapidly when M^o 
goes below ~ 200 GeV, this is because for such M^o values, the Higgs boson mass is just above 
the production threshold of the Z boson bosons, hence the result is strongly suppressed by the 
available phase space. As in the decay of — )> WW, for tan/3 = 5 the fermionic and sfermionic 
contribution comprises the dominant part of the 0{a) corrections, while for tan/3 = 30 the 
contribution from other sectors becomes important. The relative difference between the effective 
one-loop results and the 0{a) ones is less than 50% in the plotted region for tan/3 = 5. For 
tan /3 = 30 the effective one- loop results including the fermionic and sfermionic loop contribution 
differ significantly from the corresponding 0{a) ones at large M^o values, but the difference 
between the complete effective one-loop result and the complete 0{a) result is smaller, with a 
relative size less than 50% when M^o < 550 GeV (it reaches 100% for M^o ~ 800 GeV). Fig. [2] 
and [3] show that although the effective Born decay width can be rather small for certain M^o 
values (especially for large tan/3), the one- loop corrected widths for different values of M^o turn 
out to be of comparable size. 

Fig. m and [S] show the lowest order and the corrected partial decay widths of — >• WW and 
ZZ as a function of tan /3 in the m^^^ scenario for M^o = 200, 500 GeV. For both M^o 
values, the partial decay width decreases with tan /3. As in previous plots, the leading contribution 
from the fermionic and sfermionic sector is from the third generation fermions and sfermions. For 
Mj^o = 200 GeV, the fermionic and sfermionic contribution comprises the dominant part of the 
0{a) corrections, while for M^o = 500 GeV the contribution from other sectors also becomes 
important. The relative difference between the (complete) effective one-loop results and the 
complete 0{a) ones is less than 60% for both M^o values in the figures. The fact that the partial 
decay width of — )• ZZ is smaller than that of — >■ WW for the same values of tan/3 and 
Mj^o is due to the presence of identical particles in the ZZ final state. 

In Fig. O we show the corrected partial decay width as well as the relative size of the radiative 
corrections for — > WW in the M^o-tan /3 plane for three different scenarios. For the size of the 
width and the relative corrections see the caption of the figure. As illustrated there, in the 
scenario the width is rather small for large tan /3 in a wide range of M^o values. It increases 
when tan/3 decreases, the relative size of the loop corrections increases rapidly with M^o and 
exceeds the effective tree-level result when M^o > 310 ~ 460 GeV depending on the values of 
tan/3. For large values of M^o (tan/3 < 30), the relative size of the loop corrections decreases 
with Mj^o and becomes negative when M^o ~ 650 GeV. In the no-mixing scenario, the corrected 
width also increases when tan/3 decreases, and the fermionic and sfermionic contributions are 
negative over a large fraction of the scanned parameter space. In the small-agg- scenario, the 
corrected decay width increases when tan /3 decreases unless when both M^o and tan /3 are large, 
where the partial decay width is significantly increased by the Higgs propagator corrections. The 
relative correction is negative in the scanned M^o-tan /3 plane except in the upper-left corner. At 
large tan /3 values, the relative size of the loop corrections increases with tan /3 and exceeds 100% 
quite rapidly. Fig. [7] illustrates the results for — >■ ZZ in three different scenarios. The results 
shown in these figures exhibit similar features to those shown in the plots for — )• WW, but 
the corresponding width is smaller due to the presence of identical particles in the final state. 
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5 Conclusions 



We have computed the electroweak 0{a) radiative corrections to the decay of the heavy CP-even 
MSSM Higgs boson to weak gauge bosons. Due to the suppression of the tree-level H^WW/H^ZZ 
coupling, the electroweak contributions to the partial decay width are significant, they can easily 
exceed the tree-level result in certain parameter space. Although the effective Born decay width 
can be rather small for certain M^o values (especially for large tan/3), the corrected partial widths 
for different values of M^o are of comparable size. We also presented the effective one-loop results 
for the partial decay widths, which include corrections beyond 0{a). The numeric impact of such 
corrections is significant for large values of tan/3 and M^o, while it is less important for small 
values of tan/3, and also for large tan/3 with small to moderate M^o values. 
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Figure 2: The lowest order and the corrected partial decay widths for H — )• WW as a function of 
M^o in the m^^-^ scenario for tan /3 = 5, 30, where "effec. Born" denotes the effective Born result, 
"3rd (all) ///" denote the results including contributions from the third (all) generation fermions 
and sfermions, "compl. 0{a)" means the result including the complete 0{a) contribution, "impr." 
denotes the effective one-loop result. 




Figure 3: The lowest order and the corrected partial decay widths for — )• ZZ as a function of 
Mj^o in the 171^^-^ scenario for tan /3 = 5, 30. 
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Figure 4: The lowest order and the corrected partial decay widths for WW as a function 

of tan P in the m^^^ scenario for M^o = 200, 500 GeV. 




Figure 5: The lowest order and the corrected partial decay widths for H — )• ZZ as a function of 
tan /3 in the m^^^ scenario for M^o = 200, 500 GeV. 
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Figure 6: Results for the decay width of WW in three scenarios. The upper row shows the 

corrected decay width (effective one-loop result). The purple region corresponds to T fjo < IMeV, 
the blue region to IMeV < Tj^o < 5MeV, the green region to 5MeV < Tj^o < lOMeV, and 
the yellow region to lOMeV < Tjjo < 50MeV. The lower row shows the corresponding relative 
correction 5 (divided by the effective Born result). The purple region corresponds to < < 5%, 
the blue region to 5% < 6 < 25%, the green region to 25% < S < 50%, the yellow region to 
50% < 6 < 100%, and the orange region to 6 > 100%, the black region corresponds to negative 
relative correction. 
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Figure 7: Results for the decay width of — ?• ZZ in three scenarios. The upper row shows the 
corrected decay width (effective one-loop results). The purple region corresponds to Fj^^o < IMeV, 
the blue region to IMeV < Fj^^o < 5MeV, the green region to 5MeV < Fj^^o < lOMeV, and the 
yellow region to lOMeV < Fj^^o < 50MeV. The lower row shows the corresponding relative 
correction b (divided by the effective Born result). The purple region corresponds to < < 5%, 
the blue region to 5% < b < 25%, the green region to 25% < 8 < 50%, the yellow region to 
50% < (5 < 100%, and the orange region to 5 > 100%, the black region corresponds to negative 
relative correction. 
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